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KINETIC  APPROACH  TO  RADIATIVE  NONEQUILIBRIUM  FLOW 
WITH  APPLICATION  TO  GAS  FLOW  LASERS 

Gao  Zhi,  Institute  of  Mechanics,  Chinese  Academy  of  Sciences 
ABSTRACT 

A  kinetic  approach  to  nonequilibrium  flow  of  lasing  gas  is 

presented.  The  author  introduces  a  new  gain  related  to  molecular 

speed  (GMS)  and  develops  an  approximate  method  of  solution. 

These  treatments  make  it  possible  to  exactly  describe  the 

interaction  between  radiative  field,  macroscopic  flow  and 

microscopic  molecular  motion.  In  the  case  of  C02  gas  flow 

lasers,  the  zero-order  approximation  solutions  of  this  theory  are 

already  satisfactory  in  that  they  are  valid  for  the  whole 

pressure  range.  The  results  of  the  zero-order  solutions  agree 

well  with  numerical  results,  and  are  in  accordance  with  those  of 

the  currently  accepted  rate-equation  theory  (RET)  in  the  high 

pressure  range.  For  zero  flow  speed,  this  theory  leads  to  the 

well-known  theory  of  non-flow  gas  lasers  [11].  One  of  the 

present  conclusions  is  specially  worth  noting,  i.e.,  when  low- 

pressure  broadening  constant  \<  0.2,  the  rate-equation  theory, 

although  the  line  shape  factor  of  the  revised  pressure  effect  was 

introduced  [4,5],  cannot  correctly  account  for  the  effects  of 

inhomogeneous  broadening.  For  example,  0.02,  Ir/Ir  are 

about  8  “'*><?»»  £=0  and  20  when  £-=1.0,  where 

as  - 

4^ the  frequency  shift  parameter,  IR  and  Ir  are  the 

dimensionless  radiative  intensities  of  RET  and  this  theory, 

respectively. 
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NOTATION 


c  is  the  speed  of  light 

Cp  is  specific  heat  at  constant  pressure 

Fi'  Fi  speed  distribution  function  of  i-th  energy -level 

particle,  and  distribution  function  of  its  equilibrium 
speed 

/.»  /?/.  -  j J”  F,dV<dV„ 

!'  -  \  X-J' dV ^ ' 

fy  distribution  function  of  photon 

G  gain  coefficient 

G>j>  gain  correlated  to  molecular  speed 

h  Planck  constant,  or  static  entropy  of  gas  flow 

J,  Jt  radiation  intensity,  saturated  strength  and 

penetrating  radiation  intensity 

*T/  kr  velocity  of  elastic  collision,  characteristic  velocity 
of  radiation 

kij,  Kj^j  velocity  of  inelastic  collision 

Li(i  =  1,  2,  3)  length  of  optical  cavity  along  the  direction 

of  the  coordinate  axis 

1,  lx  ly,  12  direction  vector  of  light  propagation,  and 

three  direction  cosines 

m  molecular  weight 

n^  particle  number  density  at  i-th  level 

p  gas  pressure 

Rj_  =  1  -  a^  -  tj_  reflective  index  a^  of  mirror  is  absorption 

rate,  while  t^  is  the  penetrating  rate 
T,  u  gas  flow  temperature  and  flow  velocity 

V.  Vrp  particle  velocity  vector  and  thermal  velocity  vector 

r.  -  [’  J"  r,4vjiv.  pumping  velocity 
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v,  ^1,1-1  light  frequency,  transition  frequency  from  i-th  to 
(i-l)-th  energy  level 

Ai>D'  AJfy  whole  widths  at  half -peak  value  for  inhomogeneous 
and  homogeneous  broadening  type  lines 
£  frequency  shift  parameter  or  transformation  coordinate 

H  broadening  parameter 

,  JL2  intrinsic  values 
P  gas  flow  density 

6  constant 

0  superscript  0  denotes  motion  along  the  gas  flow 

direction,  the  initial  position  of  laser  oscillation 

I.  Introduction 

In  the  study  of  radiation  in  equilibrium  flow  on  the 
interaction  between  radiation  and  gas  flow,  emphasis  is  placed  on 
the  particle  characteristics  of  radiation,  but  there  is  no 
consideration  of  wave  motion  structure  [1,4].  Generally,  the 
study  can  be  divided  into  two  categories:  1)  Kcn  >  Kr,  this  is 
the  situation  for  the  study  of  the  physical  gas  dynamics  [1-3]; 
Kcn  and  Kr  are,  respectively,  the  characteristic  velocity  of 
intermolecular  elastic  collision  and  radiation  transfer.  In  this 
situation,  the  molecular  distribution  of  quantum  energy  levels  is 
controlled  by  the  collision  process.  The  radiative  transfer  of 
the  energy  is  a  nonequivalent  process.  2)  Kcn<Kr,  this  is  the 
situation  with  the  existence  of  anti-Boltzmann  distribution  and 
laser  emission  as  the  fe  ;ure,  such  as  gas  flow  lasers  [4,9].  To 
calculate  the  motion  of  laser  medium  gas  and  its  radiation 
properties,  generally  the  sets  of  simultaneous  equations  of  fluid 
dynamics,  radiation  transfer,  and  velocity  equations  are  solved 
simultaneously  [4-6]  (the  set  of  velocity  equations  describes  the 
variation  of  the  Boltzmann  constant  of  the  energy  level).  For 
convenience,  this  is  called  the  rate  equation  theory  (RET).  In 
RET,  it  is  assumed  that  particles  at  different  velocities  at  the 
same  energy  level  can  react  with  a  monochromatic  radiation  field. 
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therefore,  the  inhomogeneous  broadening  effect  cannot  be 
correctly  reflected.  As  pointed  out  by  authors  in  reference  [6], 
RET  is  not  suitable  to  be  used  in  the  situation  of  medium  to  low 
gas  pressure,  but  is  suitable  only  for  cases  of  high  gas 
pressure. 

From  radiation  theory,  we  know  that  only  frequency-resonant 
molecules  [7]  (the  doppler  frequency  for  absorbing  or  induced 
emission  molecules  approaches  the  frequency  of  the  radiation 
field)  can  have  direct  interaction  with  a  monochromatic  radiation 
field.  Therefore,  if  the  monochromatic  radiation  field  is  very 
intense  and  the  gas  pressure  is  low,  that  is,  the  homogeneous 
broadening  is  predominant,  the  frequency-resonant  molecules 
within  an  energy  level  will  be  surplus  (absorption  situation)  or 
insufficient  (emission  situation).  That  is,  the  velocity 
distribution  function  of  the  energy  level  will  be  protruding,  or 
burning  a  hole  [7].  It  is  not  possible  in  RET  to  separate  the 
molecules  between  the  frequency-resonant  molecules  and  those 
molecules  unable  to  directly  affect  the  radiation  field  because 
of  excessive  doppler  frequency  shift.  It  is  necessary  to  explore 
the  more  rational  model.  The  Lamb  theory  [8]  and  its  extension 
make  possible  an  ideal  treatment  of  the  inhomogeneous  broadening 
effect  in  a  situation  in  which  gas  properties  do  not  vary  with 
time  and  space.  This  article  explores  some  aspects  of  kinematics 
of  gas  flow  lasers  in  the  situation  when  gas  properties  vary  with 
the  flow  direction  distance.  The  kinematic  equations  describe 
the  variation  of  rate  distribution  function  of  energy  level 
particles.  In  a  study  of  kinematics,  the  interaction  among  gas 
particles  in  which  the  radiation  field  and  the  thermal  and 
macroscopic  motion  can  be  well  described.  However,  it  is  very 
difficult  to  solve  the  set  of  simultaneous  equations  relating  to 
coupling  between  the  flow  field  and  the  relaxation  process,  on 
the  one  hand,  and  radiative  transfer,  on  the  other.  In  this 
article,  a  new  physical  concept  is  introduced,  concerning  gain 
relating  to  thermal  molecular  velocity;  moreover,  an  approximate 


4 


solution  method  is  developed  to  overcome  the  difficulties.  This 
concept  is  quite  effective  for  use  in  gas  flow  lasers. 

II.  Model  of  Kinematics 


(2.1) 

(2.2) 

(2.3) 


1.  Fundamental  set  of  equations:  regarding  the  quantum 
energy  levels,  the  set  of  kinematic  equations  and  the  steady- 
state  radiative  transfer  equation  are,  respectively,  as  follows: 

PIl  _  (V  +  W,)gr adF,  -  r,  +  kAF*  -  F.)  +  (  -  V^iV^iY* 

dt  J 

-  j  Kj„_,F1WVr1_,  -  J  K^F.dV,^  +  j  K,-i.,Fl-i8(V’  -  V,)dV.-idV’ 

+  |0  |0  ~  t*iF,)dvdQ 

—  Bj.i-iCt,  F ,  —  B -i F i-^)dvdQ 

r/grad/,  =  /,  |  d>,.i-i(B,„-la,Fi  —  B,-t.,a,-|F i-t)dv 

-  if  {[*  - (i  +  7  w  • ')]' + (^) } 

The  set  of  equations  (2.1)  describes  the  relaxation  process  of 
the  initial  inequilibrium  distribution  toward  the  local- 
equilibrium  Boltzmann-Maxwell  distribution.  In  (2.1),  the 
elastic  collision  integration  was  replaced  by  the  B-G-K  model 
[1].  The  inelastic  collision  term  is  expressed 

phenomenologically.  In  the  inelastic  collision  and  the  radiation 
terms,  only  a  mono-quantum  jump  is  considered;  generally,  a 
multiquantum  jump  can  be  neglected.  Radiation  pressure, 
spontaneous  radiation  and  the  contribution  made  by  scattering  are 
also  neglected.  Ki+i,i  indicates  that  the  i-th  energy  level 
raises  a  particle  to  the  given  velocity  category;  the  (i+l)-th 
energy  level  simultaneously  loses  the  collision  transfer  velocity 
constant  of  a  particle;  indicates  that  the  (i-l)-th  energy 

level  increases  a  particle;  the  i-th  energy  level  given  velocity 
category  simultaneously  loses  the  transfer  velocity  constant  of  a 
particle.  and  are  related;  this  relationship  can 

be  derived  from  the  principle  of  detailed  balance;  F^  is  the 
pumping  term,  such  as  electron  excitation,  photoexcitation  and 
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excitation  by  chemical  reaction ,  and  so  on. 


2.  Gain  relating  to  molecular  velocity  ( GMS ) .  When  the 

characteristic  radiation  velocity  K,(Kr  is  greater 

than  the  characteristic  inelastic  collision  velocity 

"»  j  and  is  comparable  with  the  elastic 

collision  velocity  kip,  in  the  energy-level  spectral  line  shape, 
only  frequency-resonant  molecules  (the  molecules  with  consistent 
doppler  frequency  of  absorption  or  induced  emission,  and  the 
frequency  of  the  monochromatic  radiation  field)  can  have  a  direct 
function  with  the  monochromatic  radiation  field.  In  the  spectral 
lines,  the  doppler  frequency  shift  of  other  particles  is 
over large,  so  it  is  unable  to  be  directly  related  to  the 
radiation  field.  Therefore,  for  the  local  deformation  of  energy- 
level  spectral  lines,  the  radiative  transfer  can  be  in 
competition  with  elastic  collision  transfer;  the  energy-level 
distribution  function  can  possibly  have  a  protruding  or  burning  a 
hole  in  local  places  [7].  To  describe  this  physical  process,  we 
introduce  the  gain  Gt,(GMS),  relating  to  molecular  velocity. 

The  definition  of  GT ,  is  as  follows: 

Gt,  —  — - —  ,  —  B, F  ,■-»)  (2.4) 

c  *Avp 

CtA.,-,  indicates  the  gain  coefficient  of  the  gap  for  unit 
molecular  velocity  and  the  unit  three-dimensional  angle. 

Integrate  the  CTi  versus  the  apparent  frequency  v'  of  the 
molecule  from  -infinity  to  infinity  to  obtain  the  homogeneous 
broadening  gain  coefficient  G^  in  the  conventional  sense 

G*.,  -  J  G^iv  (2.5) 

As  is  the  case  in  gas  kinematics,  approximating  Gh  can  simplify 
the  problem  [1,2],  approximating  Gr,  can  possibly  simplify  the 
problems  relating  to  thermal  molecular  motion. 

3.  Approximate  solution  method.  GT,  is  the  function  of  the 
molecular  apparent  frequency  v' ,  therefore  in  Eq.  (2.1),  G Tl  can 
be  removed  outside  the  signs  of  the  double  integral  of  frequency 
v  and  the  solid  angle  fi  ;  if  the  GTl  is  considered  as  a  function 
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of  and  (Ff— F.)  ,  then  from  Eq.  (2.1)  the  approximate 

solution  of  Gr,  can  be  obtained.  On  the  other  hand,  we  solve  the 
following  double-parameter  perturbation  solution  of  Eq.  (2.1): 


f.- 2  £(—)'( fr)V 

/•D  *-»  Kt' 


(2.6) 


Obviously,  F?  is  the  Maxwell  distribution  function;  ¥f  is  the 
Chapman-Enskog  solution.  By  substituting  the  approximate 
solution  of  G T,  and  the  perturbation  solution  Fj_  in  Eq.  (2.2), 
the  sc.ation  of  the  radiative  transfer  equation  (2.2)  can  be 
obtained.  By  using  the  solution  of  Eq.  (2.2),  a  simultaneous 
solution  of  the  equation  of  macroscopic  motion  (the  moment  of  the 
kinematics  equation),  the  flow  field  variate  p,  T  and  the  flow 
velocity  can  be  obtained. 


By  utilizing  the  above-mentioned  concept  and  methods,  the 
following  problem  can  be  handled:  1)  the  situation  of  weak 
radiation,  2)  the  situation  in  which  the  time  and  space  variation 
of  the  velocity  distribution  function  is  secondary,  and  3)  the 
situation  with  high  radiation  intensity  and  discrete  frequency 
with  a  finite  number  of  discretenesses.  For  the  situation  of  C02 
gas  flow  lasers,  the  zero-level  solution  is  better  than  the 
results  of  conventional  rate  equation  theory  (RET)  [4-6. 


III.  Gas  Flow  Lasers 


In  CO 2  gas  flow  lasers,  the  light  beam  direction  is 
perpendicular  to  the  flow  direction  (refer  to  Fig.  1);  the  flow 
in  the  optical  cavity  is  approximately  one-dimensional  in  nature. 
The  effect  of  viscosity  can  be  neglected  and  the  pumping  function 
is  uniform  and  continuous.  The  molecular  relaxation  of  a  C02  gas 
mixture  is  consistent  with  reference  [5].  The  relaxation  model 
is  composed  of  five  energy-level  croups  (refer  to  Fig.  2); 
therefore,  five  velocity  distribution  functions  are  required. 

That  is,  F,(i  —  o,  1,  2,  3)  and  F0',  0  indicates  the  ground  state 
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of  CO?  vibrations;  1  and  2  indicate,  respectively,  the  CO2 
symmetric-bending  and  nonsymmetric  vibrational  models.  0'  and  3 
are  the  ground  state  and  oscillation  model  of  the  diatomic 
molecule.  The  P-branch  laser  jump  of  C02  occurs  between  the 
vibrational-rotational  energy  level  (0,  0,  1;  j)  and  (1,  0,  0; 
j+1);  j  is  the  number  of  rotational  quantum.  Overlapping  each 
vibration  energy  level,  a  series  of  rotational  energy  levels  are 
not  shown  in  Fig.  3;  however,  the  effect  of  rotational  energy 
levels  has  been  absorbed  in  factors  and  a2. 
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Fig.  1.  Optical  cavity  and 
coordinate  system 
KEY;  1  -  direction  of  gas 
flow  2  -  mirror  of  optical 
cavity  3  -  laser  output 


Fig.  2.  Relaxation  model 
of  molecular  system 
KEY:  a  -  energy  quantity 
b  -  C02  ground  state 
c  -  diatomic  molecule 


The  laser  beam  is  parallel  to  the  y-axis,  l’VT~VTf  .  The 
x  and  z  components  cf  the  thermal  molecular  motion  do  not  affect 
the  radiation  field;  therefore,  we  can  integrate  Eq.  (2.1)  with 
respect  to  vT,  and  VT,  to  obtain  the  kinematics  equation  of 


“  ^  “  r>  +  Klifi  —  I >)  +  Kuft  —  {fit I  +  /,<#>«( BuaJ,  —  B.jO,/,) 

“  ^  “  7>  +  Krdi  —  ft )  +  {itfi  —  {uit  —  {lift  —  f4>iXBuaift  — 

-  ^  -  r>  +  U/I  -  10  -  Wi  +  Ut 

(/.  -*■/,  +  l,)JVTy  -  const  (/„•  +  U)dv t,  -const 

of- j  \‘_m  ♦«(  W,  -  B^lOdv 


(3.x) 

(3.2) 

(3.3) 

(3.4) 


(3.5) 

In  deriving  Eqs.  (3.1) -(3. 3),  the  following  approximations  were 
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adopted : 


{_„  L  V^Ti6F'dV’‘dV‘>  “  (((  VTtlF,dV„dV., ) 


+  — 
dy 

*L 
”  ay 


(II  V*«F*V"*V-)  +  £(JJ VT',F.dV„dVt) 


j  K.^.F.+MV’  -  v,)dV^d\r  «,  4,+1.,F)+1 
I  -  V,)dV,-ldV’  **  {,-t.iF,-,  «  j  K.^F.dV^ 


(3.6) 

(3.7) 

(3.8) 


It  was  expressed  in  Eq.  (3.6)  that  the  variation  of  along  the 
y-axis  was  neglected.  Eq.  (3.7)  is  another  phenomenological 
expression  of  the  inelastic  collision  term;  this  actually  is 
consistent  with  the  phenomenological  expression  in  the  last 
section.  As  expressed  in  Eq.  (3.8),  the  transition  velocity  from 
the  i-th  to  the  (i-l)-th  energy  level  is  greater  than  the  reverse 
process,  that  is,  the  transition  velocity  from  the  (i-l)-th  to 
the  i-th  energy  level.  However,  \uf,  can  be  comparable  to 
4»F,  .  This  is  because  the  transfer  between  the  diatomic 

molecule  vibrational  mode  and  CO2  asymmetric  vibrational  mode  is 
near-resonant . 


On  a  mirror  surface,  radiation  satisfies  the  following 
boundary  conditions: 

y  -  0,  -  R,/,-;  y  -  L„  JZ,  -  RiJt,  (3.9) 

In  the  equation,  J+  and  J-  are,  respectively,  the  radiation 
intensities  of  positive-  and  negative-direction  propagation  along 
the  y-axis:  /  -  chvf.. 


IV.  Solution  Procedure 
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In  the  case  of  monochromatic  radiation,  for  computational 
convenient  Gtj.  is  rewritten  as 

Gr  —  ~  **(  BumJt  —  (4.1) 

By  utilizing  Eqs.  (3.1),  (3.2),  and  (4.2),  we  can  derive  the  fact 
that  the  control  equation  of  G<p  is: 

'i”1  +  +  hku  +  ■+■  ^-jcr 

i 

-  txr ,  -  ttr,  +  ^  'Mfi  ~  /.)  +  *iknt i  ~  'i(4«  +  ‘An  ~  hk»)f*  .  M  _  . 

i-i  (4.2) 

in  the  equation,  |  —  f  J.  jx 


(*  ”  f>  /i>  /j  “*  *i/>  +  *o6t>  / 1  “  /,/,  —  i®<7, 

<0  “  f  l  (  B  ji«,  +  t,  Bixa,(  &2iai  +  B|,a,)_l ,  J 

(4.3) 

According  to  Eq.  (2.4),  the  dual  parameters  of  fj^  can  be  expanded 
into 


/•- s(rV)',J 
>  -0  X 


(4.4) 

By  substituting  Eq.  (4.4)  into  Eqs.  (3.1)-(3.3)  and  into 

Eq.  (4.1),  we  obtain  the  result  that  the  zero-level  solutions  of 


f^  and  Gt  are 


/r  —  /?  — c/r  —  /?)«*  ~-u  ( ~ o 

G?  -  ^-(Bnajf  -  **  (B„aS,  - 

c  c 

f^  is  the  Maxwell  distribution,  that  is, 

/*  f>  -  „,M(r)  («  -  i,  2,  3) 


(4.5) 

(4.6) 


(4.7) 


By  integrating  the  molecular  apparent  frequency  F'  with  respect 
to  Gt° ,  we  obtain  the  generally  adopted  gain  coefficient  [4,5]. 
From  Eq.  (4.2),  we  obtain  a  more  precise  approximate  solution  (^j 
for  Gt  than  G^&  .  G*£  is  called  the  semi-order  solution.  Thus, 

first  we  obtain  n^,  that  is,  the  specifically  expressed  equation 
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of  f^.  By  processing  Eqs.  (3.1)  and  (3.2),  we  obtain 

-®Ji  +  a ,  ^  +  JJi  —  '«(|*  —  In)  +  AiG%  +  A* 
dp  d| 

°*  *  *  (4.8) 

In  the  equations,  A^  and  are  functions  of  k^j  and  Yj_;  k^j  is 
also  a  function  of  p  and  T.  As  expressed  in  experimentation  and 
analysis,  then  ^oc^T'  (0  <  fi  <  •  Therefore,  solving  for 

Eq.  (4.8)  should  be  done  simultaneously  with  the  set  of  gas 
macroscopic  motion  equations.  To  obtain  the  approximate  solution 
of  Eq.  (4.8),  the  following  mathematical  transformation  is 
introduced: 


;  —  dx'  M  ~ ***** 

a  /-_a_  y  eV7  a 

as  ” v  oc  ’  si*’"  ac>  si  ac 

(4.9) 

By  suostituting  Eq.  (4.9)  into  Eq.  (4.8)  and  neglecting  the  small 

U 

value  term  — j=  ,  we  obtain 


i M 


JLf±  +  In  *iln  ■+■  fil««  a/?  _j_  ^  '«(!* —  In)  acr  ,  .  i  -7  . 

""  77  ac  +  n  -77—  ac  +  7*  m:  h  ' 

_a_fi  4.  I»  +  •,ila  +  hk>o  a/j  ,  ,*  *»|u  sc,-  ^ 

»{■  7;  ac  s  777* 

+  ^2><  +  P 

V  ^  In 

(4.10) 

In  Eq.  (4.10),  the  coefficient  of  the  first  order  partial 
derivative  and  the  intrinsic  values  (minus  signs)  satisfy  the 
following  relationship 
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—  /_i*_  +  f  i  4.  I'A*  *,  const 

a*  »  M»*  '  M*'  » 


*»  **» 


(4.11) 


The  intrinsic  value  (not  related  to  the  ^-approximation)  is  equal 
to  a  constant.  In  the  following,  we  discuss  the  solution  **  i}  ; 
the  solution  of  ■■  1,  can  be  extrapolated  in  our  discussion. 


The  relationship  between  G®  and  f  :  by  integrating  Eq.  (3.5) 
with  respect  to  y  and  utilizing  the  radiation  boundary  condition 
(3.9),  we  can  derive 

£  J|*  -1  l„K,R,  (4.13) 


The  reflectivity  (of  the  mirror)  (i  =  1,  2)  does  not  vary  with 
x  [4-6],  therefore  we  can  generally  assume:  in —  r*cin /??/?*, 
here  5  is  a  constant  or  is  equal  to  0.  In  the  kinematics 
equation,  the  thermal  velocity  Vr,  is  not  related  to  the  space 
coordinates.  Therefore,  finally  we  have 


G’rOClnR.R,  -  e^laRlR] 


(4.14) 


Besides,  because  of  1, 
therefore  we  have 


0(0 


_L 

i,  0C 


0 


( 


uF 

hi  s/  F 


)  «  0(F), 


f 


12 


(4.15, 


By  utilizing  Eqs.  (4.14)  and  (4.15),  integrate  Eq.  (4.12)  and 
thus  we  obtain 

r. -a,  +  »).:] 

i  .|(/4*,-i,a)  —  I4  1  L  V  f* 


Ia*-  +  ( 

'  1  _  ii ka  _ 

4 

1—*)* 

•sl  M(T) } 

^  V  r*i 

1  kb 

J  J 

ft  —  +  HU'iGt  +  X)  “IT  {,»G°f~4C[^=  (1,  +  s)u/°, | 

( -  lU+'j-l.V  */  *■>  1  lV  f1 


In  the  equation 


r,  -  r,A/ (7) 


(4.16) 


n _ 1 _  [_L  -f.  CAlO 

(l,  +  5)(ij  +  «)  1/,  y/  ft 

_  1 _ +  ka  \ 

(1,  +  «)(l,  +  S)  \y/  ft  ifa  / 


(4.17) 


By  substituting  Eqs.  (4.14)  and  (4.16)  in  Eq.  (4.2),  the 
semiorder  of  G-p 

G*r-  G..M(T)  (/+  — 2 — 

In  the  equation,  the  specific  expression  7  —  JL  ,  Js  and  Gc 
be  referred  to  (5.3) 


(4.18) 


Flow  field  solution:  by  utilizing  the  solution  f®  and  G^.,  as 
well  as  Eqs.  (4.16)  and  (4.18),  the  one-dimensional  nonadiabatic 
flow  equation  can  be  obtained;  this  is  the  solution  of  the  moment 
equation  of  the  equations  (3.1)  through  (3.4): 

puA  —  const 
pu1  +  p  —  const 

2  J«  8£  W  J—  pV  t* 
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(4.19) 


h-  C,T+  f-P^T 

t  P  m 

In  the  equation  ^  is  energy  of  the  energy  vibrational  level. 
Here  we  obtain  the  zero-level  approximate  solution  of  the  set  of 
kinematics  and  radiation  transfer  simultaneous  equation  (3.1) 
through  (3.5)  for  the  C02  gas  flow  laser.  Given  lnR-~2,  Vi  and 
the  initial  conditions,  we  can  determine  the  twelve  unknowns  P, 
u,  p,  T,  h,  J,  f°  (i  =  0,  1,  2,  3),  fg  and  G$  from  12 
relationship  equations  (3.4),  (4.3),  (4.5),  (4.12)  and  (4.18). 

In  the  following,  some  useful  relationships  are  derived. 


5.  Gain,  Intensity  and  Power 


1.  Relationship  between  gain  and  intensity:  by  integrating 
the  solution  (4.18)  of  G^  with  respect  to  v '  (apparent  frequency 
of  molecules),  we  obtain  the  relationship  between  gain  and 
intensity,  namely 

G  —  (  GUv'  —  G ••&(£>  *?»  D 
1  +  7 


Av0 

J, 


-  -p-  +  xdi?  +  g-u,,  Ji2&s.  +  j!iLw 

V  p  D>‘a> 1  Vp  Vuju 


(5.1) 


in  the  equation 

Xt.,.7)-  f  _ ^ _ it  <5’2> 

«  -  VtTi,  ,  -  i(-’  ~ ».)  t/ri,  +  ±v,) 

Avd  \  c  V 

•  AV-n-  V'lnl 


*Avn  hvs/ n  '  y  n 

-.J, . + 

fit" 7={r‘+  r‘-hr-’tr')+  S  -"—{(-rL  | 

.  r; -*,*)  + 

V  <i»  y'  ft  <-  4i»  v/  p  '  */  M 
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(5.3) 


H  —  +  *»iu  — 

when  v=vq  (that  is,  the  optical  frequency  and  the  linear  center 
frequency  are  consistent),  Eq.  (5.2)  can  be  simplified  as 

c  —  exp [ t;J(  1  +  7)]  •  [1  +  erf  («?n/i  4  7)1  ^  5 

V  i  4  I 

Eqs.  (5.2)  and  (5.4)  are  adaptable  when  gain  is  equal  to  loss. 


2.  Intensity:  According  to  the  definition  and  relationship 

(3.5),  the  penetrating  intensity  is  derived  as 

J,  -  tJ 4  -  U  \/ *,  +_h  y/Jt)Lj,  Gu?*(£,  n,  H  (5.5) 

(>/*»  4  -y/RtRa)  1  +  7 


If  at  one  end  there  is  a  totally  reflective  lens  without  any 
loss,  that  is,  R2=l,  and  on  the  other  end,  there  is  the 
penetrating  output  lens,  for  the  situation  of  mainly  v  =  vq,  as 
well  as  homogeneous  and  inhomogeneous  broadening,  we  obtain, 
respectively,  the  following: 


J, 

J, 


-^(go.L2 
«1  4  /,  \ 

-*£-(  "i'L, 

a,  4  /,  v 


(g8.L2+ 

(„,iLi^L+  J_ln  J?,\ 

'  G  2  J 


(5.6) 

(5.7) 


When  u=0,  and  p  and  T  are  constant,  the  above  formulas  are 
simplified  into  a  well-known  relationship  [11]  of  gas  (not 
flowing)  lasers;  however,  we  should  pay  attention  to  the 
distinction  between  them.  Here,  G0n  "  G0n(£) • 


3.  Power:  power  can  be  obtained  by  integrating  with 
respect  x.  For  the  situation  v=Vq  and  one  end  output,  we  obtain 
the  output  power  P  as  follows  from  Eq.  (5.5): 

=  V  *  “Pj^L1  +1)1  [  1  -  erf  (,  STTbUZ  (5*8) 

1-1  (a,  4  <,)%/  P  v  1  4  7 

In  the  equation,  VD=L1L2L3  with  mainly  homogeneous  and 
inhomogeneous  broadening,  Eq.  (5.9)  can  be  converted  into 


a,  +  t, 


yin*?) 


(5.9) 
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H  - 


(5.10) 


4,  +  I, 


In  the  equation 

n  -  f  ic,  c;.  -  ({‘  « [ \1' -^t=- «r 

CL'  _  (£‘  «Gi.Ay,»  I'1 

G*  Jo  U*  J 

PL2/Vd  is  the  penetrating  radiation  intensity  by  averaging  the 
output  length  area;  Eqs.  (5.9)  and  (5.10)  are  consistent  with  the 
corresponding  power  relationship)  [11]  of  the  gas  (not  flowing) 
laser.  However,  here  Gon  and  G*2/G*  are  the  average 
quantities  in  the  flowing  direction.  By  utilizing  Eqs.  (4.15) 
and  (5.3),  we  can  derive  the  approximate  explicit  expression 
equation  as  1%  Gon  and  G*2/G*. 


6.  Analysis  and  Discussion 


1.  Comparison  with  the  exact  numerical  solution:  refer  to 
Table  1  for  parameters  used  in  the  exact  numerical  solution;  the 
corresponding  broadening  parameter  r?  is  equal  to  2.5;  this  is  the 
situation  mainly  of  homogeneous  broadening.  The  exact  results 
are  obtained  from  the  simultaneous  solutions  of  the  fluid  dynamic 
equation  and  the  rate  equation.  In  the  calculations,  the 
condition  is  used  in  which  gain  is  equal  to  loss.  The  results  of 
the  approximate  solution  and  the  exact  numerical  solution  match 
quite  closely  (refer  to  Figs.  4  through  6).  We  should  point  out 
that  the  double  integration  item  in  the  energy  equation  solution 
(4.18)  can  be  integrated  by  the  same  method  as  power  integration. 


2.  Comparison  with  the  rate  equation  theory  (RET):  usually 
simultaneous  solutions  of  RET  [4,5]  are  obtained  for  the  set  of 
fluid  mechanics  and  the  rate  equations.  For  comparison,  in  the 
following  we  briefly  derive  the  results  of  RET  corresponding  to 
Eq .  (5.1). 
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The  set  of  rate  equations  of  the  C02  laser  gas  mixture  are: 


“  “■  r,  +  +  4u)*i  +  kvCj 

*  ”  y»  +  ~  (k*  +  ku)*t  +  —  kvCj 

“  “  r  J  —  +  4u»» 


(6.1) 


TABLE  1.  Calculation  Conditions  and  Parameters  at  Inlet 
of  Optical  Cavity 


*i  +  i,  —  1,  <■  —  O.M 

»,  _  71> 

~  NT 


(***)  i 


X  - InR.  -  5  X  1°-'  (*#-)  j 

^  a 

«,  -2.8X10-'*  b 

,.-,,-4.7X10-"  (fcft/tt?)  b 
1.9  X  10-'» 

■■  —  2.92  X  10-»  (jg)  C 

0-1.4X10’  d  (*»/$.  *K) 

2*i/N-0.5 


Nt-  7.7JX10**  (**/■*')  e 

wj  —  2. jo  x  10“  (»?/«*')  e 

VJ-  3.83  X  10” 

V  -  <».M  X  10" 

V- 4.0x10*  (aa/**1)  f 
«*-  1.4  X  10*  (M*/#)  g 
7"*  —  300  °K 

K%  =  5.68  x  10-  (*Jr')  h 
jq,  =•  2.77  X  in* 

iq.  =  1.02  x  io* 


(AT,.,  *,.)  -  (/q„  **„  Kl)  ±  y/i 

JC„-  x„«p 


jq,  =  4.02X10’ 
CO./N./Ht  =  1/4/5 


KEY:  a  -  constant  b  -  (ergs  per  particle) 
c  -  (gram)  d  -  (ergs/gram,  °K)  e  -  (particles 
per  cubic  centimeter)  f  -  (dynes  per  square 
centimeter)  g  -  (centimeters  per  second) 
h  -  (second-1)  i  -  (square  centimeter) 
j  -  (centimeter-1) 
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Fig.  3.  Variation  of  n3  and  Fig.  4.  Variation  of  pressure 
nfo  with  f  and  temperature  with  f 

Legend:  -  exact  numer-  Legend:  -  exact  numerical 

cal  solution  solution 

A  o  is  the  approximate  solu-  o  A  is  the  approximate  solu¬ 
tion  in  this  article  tion  in  this  article 

A  n3  o  njj  o  pressure  A  temperature 

KEY:  a  -  pressure  p  x  10-4  dy 
per  square  centimeter^ 
b  -  temperature  TxiO-*  °K 


1  2 
c 


Fig.  5.  Variation  of  power  density 
with  f 

Legend:  -  exact  numerical  solu¬ 

tion:  o  approximate  solution  in  this 
article 

KEY:  a  -  centimeter-3 


After  introducing  the  linear  vector  [4,5]  of  the  revised  pressure 
effect,  the  gain  coefficient  is 


G 


_ l 

*Avv 


0(f»  *lf  0)( 5aa^»j  — 


(6.2) 


By  a  derivation  that  is  similar  to  that  in  section  4 
utilizing  the  condition  that  gain  is  equal  to  loss, 
is  derived: 

q  «_  £*■$(£>  Hi  o) 

l  +  *j.  °) 


,  and  by 
the  following 

(6.3) 


This  equation  and  the  references  [5,6]  have  the  same  results, 
applicable  when  gain  is  equal  to  loss.  When  high  pressure  rj  '»! , 
Eqs.  (6.3)  and  (5.1)  are  of  the  same  order  of  magnitude.  In 
cnese  two  theories,  Gon  and  loss  G  are  the  same.  Therefore,  from 
Eqs.  (5.1)  and  (6.3),  we  derive 

j "  _  l  +  Ik _ 1  (6.4) 

4>(£»  7  k)  °) 

We  can  see  in  all  the  possible  values  of  £  and  r]  ,  the  intensity 
IR  of  RET  are  greater  than  the  intensity  in  this  theory;  refer 
to  Fig.  6;  refer  to  Fig.  7  for  further  explanations.  In  the 
figure,  by  using  I  and  £  as  parameters,  and  given  the  variation 
relationship  of  G/G0n  with  rj  ,  all  curves  in  this  theory  are 
situated  below  the  corresponding  RET  curves;  all  curves  in  the 
two  theories  are  situated  below  the  homogeneous  broadening  limit 
curves.  This  explains  the  effect  of  RET  on  low  pressure;  that 
is,  the  estimation  of  the  effect  is  insufficient  for  the 
inhomogeneous  broadening  effect.  For  the  situation  of  the 
broadening  parameter  ^<0.2,  it  is  necessary  to  adopt  the  results 
of  the  kinematics  theory. 
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(1=0) 


Fig.  7.  Variation  of  G/G0n  with 
7  ($-0.5) 

Legend:  -  this  theory 

_  _  _  _  RET  theory 

_ * _  homogeneous  broadening  limit 


20 


3.  Comparison  with  the  gas  (not  flowing)  laser  theory:  the 

0 

apparent  dependence  on  £  of  parameters  such  as  f^  is  in  the  form 
of  index  ,  such  as  Eq.  (4.16).  Therefore,  when  1>C  >  I  , 

the  relationships  between  Eq.  (4.16)  and  (5.3)  can  be  simplified 


as 


+  ‘»G% 

I+0.+ 


_ r 

,  + 


q.  ~ 


20  _  T'l  i  V  „  _l_  fo^uGr  _  (  ^  j.  ^io\ 

n  T*  ~ 4"  1  +  n+aW 

J  _  chvs2  Kllktokl*  +  L  ^31^33  +  C^n  +  ^3  +  4ji)^1o1^  V  ^ 

2  Bjia,  +  (^n  +  tik*  +  V~f* 

Gul  -hv  ( r;  +  r;  -  ^  2  y-) 

(6.5) 

In  the  equation,  the  first  order  term  of  5  is  retained,  and  the 


second  and  higher  order  terms  of  6  are  neglected;  (it  can  be 
proved  that  6«Kl).  When  the  reflectivity  of  the  mirror  does  not 
vary  with  x,  that  is,  6=0.  Eq.  (6.5)  and  the  corresponding 
equation  (5.1)  are  just  the  familiar  relationship  [11]  of  the  gas 


(not  flowing)  laser.  It  is  apparent  that  the  well-known 
relationship  [11]  of  the  gas  (not  flowing)  laser  is  a  special 
case  of  this  theory  when  u  =  0  or  2*>1 .  However,  it  should  be 
noted  that  Eq.  (6.5)  is  suitable  for  the  case  when  the  gas 
properties  vary  with  the  flow  direction.  From  x  =  0  satisfying 
the  relationship  (6.5),  the  gas  flows  past  a  distance  xp  as 


** 


lu 


(6.6) 


7.  Conclusions 

Results  of  the  approximate  theory  in  this  article  are 
applicable  to  the  entire  pressure  range;  the  approximate  results 
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and  the  results  of  exact  value  match  quite  closely.  At  high 
pressure,  the  results  are  consistent  with  the  rate  equation 
theory  that  is  generally  used.  The  familiar  relationship  [11]  of 
the  gas  (not  flowing)  laser  can  also  be  obtained  as  a  special 
case  of  the  result  of  this  article.  This  illustrates  that  the 
present  treatment  of  the  kinematics  theory,  the  introduction  of 
gain  related  to  molecular  velocity,  and  the  corresponding 
approximate  solution  method  can  serve  in  relatively  exactly 
calculating  the  macroscopic  and  microscopic  motions  of  the  gas, 
as  well  as  the  interdependent  properties  of  the  three,  including 
the  radiation  field. 

This  article  was  received  in  November  1980. 


Footnotes : 

1.  This  article  was  circulated  in  the  two  following  cases:  the 
Second  All-China  Fluid  Mechanics  and  the  First  Asia  Fluid 
Mechanics  Conference  at  bangalore,  India,  in  December  1980. 
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